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ABSTRACT 

We present the first detailed spatio-kinematical analysis and modelling of the plan- 
etary nebula HaTr 4, one of few known to contain a post-common-envelope central star 
system. Common envelope evolution is believed to play an important role in the shap- 
ing of planetary nebulae, but the exact nature of this role is yet to be understood. High 
spatial- and spectral- resolution spectroscopy of the [O ill] A5007 A nebular line ob- 
tained with VLT-UVES are presented alongside deep narrowband Ha-f [Nil] A6584A 
imagery obtained using EMMI-NTT, and together the two are used to derive the three- 
dimensional morphology of HaTr 4. The nebula is found to display an extended ovoid 
morphology with an enhanced equatorial region consistent with a toroidal waist - a 
feature believed to be typical amongst planetary nebulae with post-common-envelope 
central stars. The nebular symmetry axis is found to lie perpendicular to the orbital 
plane of the central binary, concordant with the idea that the formation and evolution 
of HaTr 4 has been strongly influenced by its central binary. 
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The current understanding of planetary nebulae (PNe) is 
that they are formed when low-to-intermediate mass stars of 
O.8M0 < M < 8M0 start to transition from an asymptotic- 
iant- branch (AGB) star to a white dwarf. iKwok et al.l 
19781 ) first put forward the 'Interacting Stellar Winds' 
(ISW) theory of nebular formation, whereby a PN is formed 
by a 'snow-plow' process. Here, a slow, dense 'superwind', 
blown while the central star is in the AGB phase, is swept 
up by a fast, tenuous wind emanating from the emerging 
white dwarf creating a compressed, thin, dense shell. This 
shell is then ionised by the hot c entral star to produce the 
visible PN. Balick & Frank (|2002l ') state that in general, the 
mass loss in these circumstances is inherently isotropic which 
should lead to the formation of a spherical nebula around 
the progenitor star. However, many PNe instead tend to 
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show more complex, axisymmetric (bipolar) morphologies 
which contradict this theory. Within the framework of the 
ISW theory, any deviations from sphericity - particularly 
bipolar shapes - requi re an asphcrical mass distribution in 
either wind phase (|Kahn fc W est 1985). This has become 
known as the 'Generalised Int eracting Stellar Winds' model 
fGISW. [Salick fc Fran3l2002l ). However, the source of this 
anisotropy is currently a matter of some controversy. 

A strong argument is starting to be built up in favour 
of bipolar nebulae being shaped by a central binary star 
syste n i (|De Marco 200 9. Miszals ki et al.l |2009| . I Jones et al] 
bOllH . lPaczvnskil l|l976i) first identified a binary pathway to 
PN formation via a common-envelope p hase, with the firs t 
such nebula discovered being Abell 63 (|Bond et al.lll97^ '). 
Here, the density contrast required by the GISW model to 
form aspherical morphol ogies arises from the ejection of the 
common envelope (CE) (|Nordhaus fc Blackman|[2006l ) . Fric- 
tional forces between the binary and CE result in the trans- 
fer of angular momentum, causing the CE to be preferen- 
tially ejected in the stellar orbital plane - the symmetry axis 
of the PN is then expected to be perpendicular to the or- 
bital plane of the binary. Therefore, observing this predicted 
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alignment in a given PN is strongly indicative that the neb- 
ular morphology has been influenced by the binary central 
star. 

Such is the strength of support for the role of binarity in 
the formation and evolution of PNe, that the working group 
PLAN-B (PLAnetary Nebula BinariefQ) has been formed in 
order to coordinate the research effort in this field. 

To date, only five PNe have been observationally 
shown to possess the predicted alignment between neb - 
ula and binary plane - Abell 63 (Mitchell et al. I l2007l). 



Abell 41 ll Jones et al.ll201(]| '). AbeU 65, (IHuckvale et al. 



Shimanskv et al.1 200917. NGC 63 37 llHillwig et al l 
Garcfa-Dfaz et al.r 2009') and Sp 1 l|jones et al.il2012l ). 
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PN HaTr 4 (PN G335. 2-03.6, a = le*" 45™ 00.2", 5 = 
—51° 12' 22.0", J2000) is known to contain a photometric 
binar y central star with a period of 1.71 days (|Bond &: Livid 
Further ^ more detailed i nvest igation by Hillwig et al. 
(in prep) and iBodman et al.l (|2012l 'l revised the period to 
1.74 days. Additionally, they find no evidence for eclipses in 
the lightcurve indicating that the inclination of the binary 
orbital plane with respect to the line of sight must be less 
than 8 5° , with the best fit to the data being from a model 
at 75° (|Bodman et al.ll2012l ). 

Based on p revious imagery - the discovery image of 
iHartl fc Tritton IT98&) and subsequent narrowband images 
of iBond fc Livid - HaTr 4 has the appearance of a 

"classical butterfly" bipolar nebula lying in the plane of 
the sky with twin lobes emanating from the central star 
in an East- West direction, as well as an indication of pos- 
sessing fainter, elongated, open-ended lobes extending in a 
North-South direction (only revealed in later, deeper im- 
agery - see section (2] figure [l]). However, imagery alone 
is insufhcient to precisely and unambiguously determine 
the nebular structure and orientation, as PNe ar e very 
prone to inclination-d ependent projection effects (|Kwokl 
l2010l : lFrank et al.lll993l '). Only by using high-resolution, spa- 
tially resolved spectroscopy is it possible to fully ascertain 
the intrinsic three-dimensional structure of a given nebula 
in order to assess whether the binary has played a role in 
the shaping of the nebula. 

In this paper, we present long-slit spectroscopy and 
deep, narrow-band imagery of HaTr 4, from which a spatio- 
kinematical model is derived, in order to investigate the re- 
lationship between HaTr 4 and its central binary star. 



2 OBSERVATIONS 

The deep HQf+[N ll] A6584 A image shown in figure[T]was ac- 
quired on 1995 April 22, using the ESQ (Euro pean Southern 
Obse rvatory) Multi-Mode Instrument (EMML lOekker et al.) 
1 19861 ') on the 3.6-m ESO New Technology Telescope (NTT) 
with an exposure time of 1800 s and seeing of 1.1". The 
image clearly shows the "bow-tie" -sh aped central nebula a s 
seen in the [Olll] A5007A imagery of lBond fc Livid l|l990ll . 
but also faint extensions to the North and South. These 
extensions are more consistent with a symmetry axis ori- 
entated roughly North-South (as opposed to the previously 
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inferred East-West axis), indicative that the "bow-tie" may, 
in fact, form the waist of the nebula. 

On 2008 July 10 and 11, data were acquired of HaTr 4 
using grating #3 on the visual-to-red arm o f the Ultravi- 
olet and Visual Echelle Spectrograph (UVES; lOekker et al.l 
12000. ') on the Keuyen Unit Telescope (UT2) of the "Very Large 
Telescope (VLT). UVES was operated in its 30" longslit 
mode with a 0.6" slitwidth (R ~ 70, 000) using a narrow- 
band filter to isolate the [O in] emission line profile and 
prevent contamination from overlapping orders. 1200 s ex- 
posures were taken at ten different slit positions (shown in 
figure [ijb)). Slits 1-5 were taken at a position angle (P. A.) 
of 7°, and slits 6-10 were taken at a P. A. of 97° in order to 
give full nebular coverage in both an East- West and North- 
South direction. The spatial scale of the observations was 

0. 17" per pixel, and the seeing was between 0.8" and 1.0" 
for all observations. 

Reduction of the spectra was carried out using stan- 
dard routines within the STARLINK software package. The 
spectra were cleaned of cosmic rays and debiased appropri- 
ately. The spectra were then wavelength calibrated against 
a ThAr emission-lamp. The data were rescaled to a linear 
velocity scale appropriate for the [O ill] emission, and cor- 
rected to heliocentric velocity, Vhci ■ The reduced spectra are 
presented in figures[2]and[3]as position- velocity (PV) arrays. 

The PV arrays presented in figure [2] are consistent with 
a nebula extended in the North-South direction, as indicated 
by the NTT Ha+[N ll] image shown in figure [U with slits 

1, 3 and 5 showing velocity ellipses. However, some data are 
missing as the slit is not quite long enough to cover the full 
North-South extent of the nebular emission (as indicated by 
the sudden cut-off at the top and bottom of the PV arrays, 
e.g. at cross-section 13.5" in figure 2(e)). The PV arrays from 
both slits 1 and 3 (figure [2j a) and (c)) appear to be closed 
at their northernmost extent (as supported by the nothern 
extension of slit 3, shown in slit 2, figure[2](b)) but open in 
the South; however, the data from the southern extension to 
slit 3 - slit 4, figure [2Id) - clearly shows that this PV array 
is in fact closed. 

The presence of two separate velocity components (cor- 
responding to the red- and blue-shifted sides of a velocity 
ellipse) is consistent with front and back walls of an ex- 
panding hollow shell extended along the slit length (i.e. in a 
North-South direction) - this orientation is perpendicular to 
the East- West orientation previously indicated by the im- 
agery presented in figure [1] The emission from the central 
(equatorial) region of this shell is significantly brighter than 
that from the North and South ends, and it is this equatorial 
region which gives the impression of an East- West 'bow-tie' 
shape, as seen in figure [ij a). 

The PV arrays presented in figure [3] are obtained from 
slit positions perpendicular to those in figure [51 They each 
show a closed velocity ellipse, confirming the presence of a 
shell extending in a North-South direction. Slit 7 appears 
to have a larger FHWM for the blue-shifted component at 
cross-section 0", whereas slit 9 appears to have a larger 
FWHM for the red-shifted component at cross-section 0". 
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Figure 1. (a) High contrast image of PN HaTr 4 taken on tlie NTT-EMMI spectrograph in H«-|-[N ll]. (b) The same NTT image 
presented at a lower contrast showing all slit positions taken on the VLT in [O III]. Slit width is 0.6". The central star is visible at 
16:45:00.2 -51:12:22.0 



SPATIO-KINEMATICAL 
RECONSTRUCTION 



astro physi c al modelling progr am 



SHAPE0 

a spatio- 



Usins the astr op 

jSteffen fc Loped [2OO6I : ISteffen et al] I2OIII ) 
kinematical model was developed in order to reconstruct 
the nebular morphology of HaTr 4 based on both the 
high-resolution Hq-|-[N 11] imagery shown in figure[TIa) and 
the high-resolution, spatially resolved [O in] spectra shown 
in figures [2] and O 

A wide range of different morphologies and orientations 
for the nebular shell and waist were applied to the model 
of HaTr 4. The nebular expansion velocity was assumed to 
be a Hubble-type flow, but the scale velocity was consid- 
ered a free parameter in the modelling ([r/ro]*k, where r is 
radius from the central star, To is the inner equatorial ra- 
dius , and k is the velocity at To). The best-fitting basic 
model, determined by eye, was an elongated ovoid nebular 
shell encorporating a thicker equatorial ring, an image of 
which is shown in figure |4] at the same scale, P. A. (7°), and 
inclination to the line of sight as the image shown in figure 
[T] and presented again here for comparison. The synthetic 
PV arrays corresponding to this best-fitting model are pre- 
sented in figure [5] The model accurately reproduces the two 
velocity components seen in all of the slit positions shown 
in figures [2] and [S] corresponding to the front (blue-shifted) 
and back (red-shifted) walls of the main nebular shell. The 
model also recreates the bright emission associated with the 
nebular waist, appearing as the two bright velocity compo- 
nents around cross-section 0" (see figure [Sfe)). 
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An ovoid model was favoured over a more cylindrical 
structure, as a cylindrical shell produced a nebula which was 
too 'straight-edged' compared to the observed imagery, and 
did not create the roughly elliptical profile followed by the 
two spectral components observed in slits 1 to 10. Similarly, 
a tighter waist was also discounted as it gave both the neb- 
ular image and spectra a 'figure-of-eight'-like appearance, 
with the nebular shell protruding from the central 'bow-tie' 
rather than it being a smooth transition i.e. it is not 'wasp- 
waisted', akin to other such close binary central star plane- 
tary nebulae (CSPNe) whose waist radi us is comparable t o 
that of its bipolar lobes such as Abell 41 (| Jones et al.ll201ol '). 
If the shell thickness parameter was set to greater than 5", 
it created broader components at the northern and southern 
tips of the model spectra than are observed in the data. A 
ring was applied to the equatorial region of the model in 
order to increase the apparent thickness of the shell at low 
latitudes. The 'bow-tie' shape could not be accurately repro- 
duced by simply enhancing the brightness of the shell in this 
region. Moreover, the increased physical thickness also leads 
to broader FWHM spectral components from this equatorial 
region due to the larger range of velocities present. 

Once the basic model was determined, its inclination 
was varied in order to assess over what range a satisfactory 
fit to the data could be found. Comparison between observed 
and model imagery constrains the inclination of the nebula 
to be between 65° and 80° - below 65° the central 'bow tie' 
becomes too extended in the North-South direction, and its 
centre becomes too narrow as less of the ring overlaps in the 
line of sight. Above 80°, the ring becomes too 'tube-like' in 
shape rather than the apparent bow tie as the model nebula 
becomes close to lying in the plane of the sky. Comparison 
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(a) Slit position 1 



(b) Slit position 2 



(c) Slit position 3 
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(d) Slit position 4 



(e) Slit position 5 



Figure 2. PV arrays showing reduced [O in] spectra from the VLT slit positions 1-5. North is up, South is down. The velocity axis 
on all plots is heliocentric velocity, D^ej. The display scales of each slit have been set individually to highlight the spatio-kinematic 
features referred to in the text. Cross-section O" defines where the central star is found, the continuum of which can be seen in slit 3. 
The continuum of a field star is visible in slit 5 at cross-section -8". 
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(a) Slit position 6 



(b) Slit position 7 



(c) Slit position 8 
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(d) Slit position 9 



(e) Slit position 10 



Figure 3. PV arrays showing reduced [O III] spectra from the VLT slit positions 6-10. West is up, East is down. The velocity axis on 
all plots is heliocentric velocity, v^el ■ The display scales of each slit have been set individually to highlight the spatio-kinematic features 
referred to in the text. Cross-section 0" defines where the central star is found. 
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(a) 75° model nebula (b) Observed image of HaTr 4 

Figure 4. SHAPE model of HaTr 4 at the best-fitting inclination of 75°, with the real image for comparison. 



Table 1. Table of the main parameters of nebular ring in the 
best-fitting shape model of HaTr 4 shown in figure |4] (errors are 
the range over which the parameter could be altered and still be 
considered a reasonable fit to the observed image and spectra). 



Parameter Value 

Semi-major axis (to outer edge) 12.65±0.5" 

Semi-minor axis (to outer edge) 7.5±0.5" 

Shell thickness 5±1" 
Equatorial expansion velocity, Vexp (r = 7.5") 13±2 km s" 

Heliocentric systemic velocity, Vsys -94±4 km s 

Inclination, i 75±5° 



between modelled and observed spectroscopy further con- 
strains the inclination to be 75° ±5°, as other inclinations 
did not replicate the observed brightness offset between red- 
and blue-shifted velocity components seen in slits 1 to 5. 

The key parameters of the best-fit model are sum- 
marised in table [1] The determined expansion velocity rela- 
tive to the central star of V"ei:p=13± 2 km s~^ in the eq uato- 
rial plane is fairly typical for a PN (|Weinbergerlll989l ) . and 
the systemic veloci ty (Kns=-94±4 km s~ ^) is in good agree- 
ment with that of iBeaulieu et aL I l| 19991 : Vsys=-97±l^ km 
s~^). The kinematical age per unit distance, as determined 
by the model parameters, is found to be 2740±420 years 
kpc~^ . Usin g his surface brightness-distance relationship, 
Frew (|2008l ) finds a distance to HaTr 4 of roughly 3 kpc; 
using this value, we derive a kinematical age of 8200±1300 
years for HaTr 4 and an equatorial radius of approximately 
0.22 pc. 

Figure [5l^e) shows the model spectrum corresponding to 
the observed PV array from central slit 3 (reproduced above 
it in figure [S^ a)) at the best-fitting inclination of 75°. The 
model PV array nicely confirms that the splitting of bright 
components shown in the spectrum is indeed due to an elon- 



gated shell possessing a thicker and brighter equatorial ring, 
as the relative vertical offset between the two bright com- 
ponents is dependent on nebular inclination. Figures [5jf), 
(g) and (h) shows the SHAPE models of slits 8, 9 and 10 
(see figure Hlb)) at the best-fitting inclination of 75° be- 
neath their observed spectral counterparts to highlight how 
the brightness variations across the velocity ellipses changes 
depending on the shell/waist overlap as a result of inclina- 
tion effects. The fact that the varations in the model do not 
exactly match the spectra highlights the asymmetry of the 
structure of HaTr 4 in reality. Indeed, in all of the spectra 
shown in figures [2] and [S] the blue-shifted emission compo- 
nent is observed to be brighter than the red-shifted with 
no obvious explanation, as the best-fitting model does not 
replicate this trend; however, possible explanations are ISM 
interaction (which may also account for the thinner nature 
of that side in figure |3jc)), or more likely that there is in- 
ternal extinction blocking emission from the back wall. 

The SHAPE- modelled nebular inclination of 75° ± 5° rel- 
ative to the line of sight is approximately perpendicular to 
the binary plane, as determined in the work of Hillwig et al. 
(in preparation) and iBodman et al.l (|2012l ) . iBodman et al.l 
(2012) state that the best-fit range for the binary falls be- 
tween 75° and 80° , but with an acceptable fit as low as 55° . 
In knowing that the PN symmetry axis is expected to lie 
perpendicular to the orbital plane of the binary and hav- 
ing seen that the inclinations of the two are in agreement, 
these results further imply that the binary is responsible for 
influencing the shaping of HaTr 4. 



4 DISCUSSION 

The imaging, spectroscopy and subsequent spatio- 
kinematical modelling presented in this paper clearly 
show HaTr 4 to have an elongated, axisymmetric mor- 
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Figure 5. Top row ((a)— (d)) reproduces the observed spectra for slit positions 3, 8, 9 and 10 (see figures [2] and [Sj ■ Bottom row ((e)— (li)) 
shows the equivalent SHAPE models, at a best fitting inclination of 75° relative to the line of sight, to show the differences between cuts 
through the major axis, equatorial ring, and the outer shell of the nebula. Slit positions are as shown in[TIb) 
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phology with an equatorial enhancement consistent with a 
nebular ring. This work demonstrates for the first time that 
the symmetry axis of HaTr 4 lies in a roughly North-South 
orientation, as opposed to the East-West orientation 
inferred from previous imagery. The "bow-tie" appearance 
of the central region of the nebula is shown to result from a 
line-of-sight inclination effect associated with the enhanced 
nebular waist, rather than a bipolar structure as previously 
believed. 

The spatio-kinematical modelling of HaTr 4 reveals 
fairly typical dynamical properties, with a kinematical age 
of order 8000 years and an equatorial expansion velocity 
of 13±2 km s~^. No evidence is found for extended emis- 
sion beyond the main nebular shell or jet-like outfiows as- 
sociated with HaTr 4, such as those found in the binary- 
centred PNe ETHOS 1 (|Miszalski et al.ll201ll'). The Nec klace 
llCorradi et all 1201 ll). Abell 63 jMitchell et all l2007h and 
NGC 6337 l|Garcfa-Di'az eTZI 120091 '). Similarly, even once 
the effects of inclination have been taken into account, no 
obvious analogue for HaTr 4 can be found amongst the other 
binary-centred PNe which have been the subject of detailed 
spatio-kinematical study. While the morphology of HaTr 4 
is axisymmetric it is not bilobed as those PNe mentioned 
above are, displaying more of an elliptical overall structure. 
Perhaps the closest m orphologicall y to HaTr 4 are A bell 63 
l|Mitchell et al.ll2007l ) and Abell 41 (| Jones et al.ll2010[ ): how- 
ever, both show marked differences to HaTr 4: Abell 63 has 
no ring-like equatorial feature, and also displays extended 
polar outflows which are absent in HaTr 4. Abell 41 does 
display an equatorial ring, but a far less extended version 
than the one seen in HaTr 4. Similarly, Abell 41 has a dis- 
tinct bilobed structure not seen in HaTr 4. This lack of direct 
analogue highlights the morphologically disparate nature of 
the current samp le of binary-centr ed PNe that have been 
studied in deta il ( Jones et al.ll201ll ). The imaging study of 
iMiszalski erahl l|2009l ) did, however, indicate that the most 
common morphological feature amongst the known sample 
is some form of equatorial ring - a feature we have shown 
to clearly be present in HaTr 4, and is probably associated 
with the ejection of a common envelope. 

The inclination of the nebular symmetry axis is found 
to be ~75°, consistent with the non-eclipsing nature of 
the binary central star. Further investigation into the cen- 
tral star system by Hillwig et al. (in preparation) and 
iBodman et all (|2012l '). indicates that the inclination of the 
central binary plane is very close to the nebular inclina- 
tion, as predic ted by theories of binary-induced PN shaping 
ijNordhaus fc Blackman 20061 . HaTr 4 is one of only six PNe 
to have had this alignment between nebular symmetry axis 
and binary plane observationally shown, a finding entirely 
consistent with HaTr 4 having been shaped by its central 
binary star. 
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